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Abstract 

The B*Btt and D*Dtt couplings have previously been derived from a QCD light- 
cone sum rule in leading order. Here, we describe the calculation of the 0(a s ) 
correction to the twist 2 term of this sum rule. The result is used for a first next- 
to-leading order analysis. We obtain gs* Bit = 22 ± 7 and Qb*d-k = 10-5 ± 3, where 
the error indicates the remaining theoretical uncertainty. 
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1. The hadronic B*Btt coupling is denned by the on-shell matrix element 

(B*°{p) n-(q) | B-(p + q)) = -g B * B ^{q-e) , (1) 

where the meson four-momenta are given in brackets and e M is the polarization vector of 
the B*. An analogous definition holds for the D*D% coupling. These couplings play an 
important role in B and D physics. For example, they determine the magnitude of the 
weak B — > it and D — > 7r form factors at zero pion recoil, i.e., at momentum transfer 
squared close to {m B — m^) 2 and {m D — m^) 2 , respectively. Moreover, the coupling 
constant Qd*d-k is directly related to the decay width of D* — > Dn. The decay B* — > Bn 
is kinematically forbidden. 

Theoretically, the B*Bn and D*Dtt couplings have been studied using a variety of 
methods^. Among these, QCD light-cone sum rules (LCSR) have proved particularly 
powerful. The first LCSR calculation of g B * BlT and Qd*d-k including perturbative QCD 
effects in leading order (LO) was reported in JIJ] . It is well known, however, that in heavy- 
light systems next-to-leading order (NLO) effects can be essential. An important example 
is provided by the heavy-to-light form factors 0, |3|, f|, ||. This calls for a NLO analysis 
of Qb* Bit and g D *Dw 

In this letter, we describe the calculation of the 0(a s ) correction to the leading twist 2 
term of the relevant LCSR. The result is then used for a first NLO evaluation. We present 
numerical predictions for gs*Bn and gr>*D-K together with an estimate of the remaining 
theoretical uncertainties. 

2. For definiteness, we focus on the derivation of the B*Bn coupling. For the D*Dn 
coupling the procedure is completely analogous. The LCSR for Qb*Btk is derived from the 
vacuum-to-pion correlation function of the vector and pseudoscalar currents involving the 
6-quark and light quark fields: 

Fp(p,q) = i J d A xe^{-K-{q) \ T{d(x)^b(x),m b b(0)i l5 u(0)} | 0) 

= F(p 2 , (p + q) 2 ) qtl + F(p 2 , (p + q) 2 ) Pll . (2) 

Inserting in the above matrix element complete sets of intermediate hadronic states car- 
rying B and B* quantum numbers, respectively, and using the definition ([1]) together 
with 

m b {B~ | bij 5 u | 0) = m 2 B f B (3) 

and 

(0 | d^b\ B*°) =m B *f B ^ , (4) 

one obtains double dispersion relations for the invariant functions F and F. In the 
following, we only need 

F(r? ( n + a) 2 ) = m2 B m B*fBfB*9B*Bn ff p h (s U S 2 )ds 1 ds 2 . * 

{P ,{P~r q) ) (p2_ m ^ )((p + ?)2 _ m | ) +yy {si _ p 2 ){s2 _ {p + q) 2y I ) 



L An overview is given, e.g., in Tab. 1 of ref. jy. 
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The first term on the r.h.s. is the ground-state contribution involving the B and B* 
masses, decay constants, and the desired B*Btt coupling. The second term results from 
the higher resonances and continuum states in the B* and B channels described by the 
hadronic spectral density p h (si,s 2 ). The integration boundary in the (si, s 2 )-plane is 
denoted by E. In general, the dispersion relation (|5]) requires subtractions in order to 
render it finite. However, these subtraction terms being polynomials in p 2 and/or (p + q) 2 
will be removed by Borel transformation^, and are therefore not made explicit here. 

In QCD, the invariant amplitude F(p 2 , {p + q) 2 ) can be calculated at p 2 , (p + q) 2 <C ml 
by expanding the time-ordered product of currents in the correlation function near 
the light-cone, that is around x 2 = (for reviews, see e.g. || 0). The result is given by 
convolutions of hard scattering amplitudes and pion distribution amplitudes. The latter 
parameterize the long-distance effects and are classified by twist. This paper mainly deals 
with the leading twist 2 term 

i 

F^(p 2 ,(p + q) 2 ) = -f 7T [ du^(u,p)T(p 2 ,(p + q) 2 ,u,p) , (6) 



o 

T being the hard amplitude, <p n the twist 2 distribution amplitude, and /i the factorization 
scale. In analogy to (|5|), we write (|) in the form of a double dispersion relation: 

J J (s x -p 2 ){s 2 - (p + q) 2 ) 
where the corresponding spectral density p is given by 

p(si,s 2 ) = ^lm Sl lm S2 F {2) (s 1 ,s 2 ) . (8) 

After equating fl5|) and ([F]) and applying Borel transformation one obtains the basic sum 
rule 

mlm B *f B f B *g B * B7r e W (-^ - ^ + JJ P h (s u s 2 ) exp (-j^ - j^j d Sl ds 2 



= JJp( Sl , 82) exp (-j^L - ^ j dsi ds 2 (9) 

with M 2 and M 2 being the Borel parameters associated with p 2 and (p + q) 2 , respectively. 

To proceed, one has to subtract from (H) the unknown contribution of the excited B 
and B* states and of the continuum states. This can be achieved approximately by using 
quark-hadron duality. To this end, the integral over p h is replaced by a corresponding 
integral over p with a suitable integration boundary E. After subtraction, one has 

r r 1 f f / \ ( m R* ~~ S l m R ~ S 2 \ , , / n 

hh * 9B ^ = ^^A\ p ^ S2) exp [-^f- + ~w~ ) 12 • ( } 

2 For the definition of the Borel transformation see, e.g., ref. 
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It should be noted that this LCSR actually determines the product of the B and B* decay 
constants and the B*Btt coupling. This will play an important role in the numerical 
analysis. 

3. The perturbative hard scattering amplitude T introduced in (H) can be expanded 
in a power series in a s : 

T(p\ (p + q)\ u, /*) = T (p 2 , (p + q)\ u) + ^Up 2 , (p + qf, u, /*) + 0(a s 2 ) . (11) 
The lowest-order term is given by [|TJ 

2 

T (p 2 ,(p + g) 2 , M ) = - 2 m \ . V2 , (12) 

ml — p z (l — u) — (p + g) 2 u 

where m;, is defined to be pole mass of the b quark. The 0(a s ) term 7\ has recently been 
calculated [0, |3| in the context of the LCSR determination of the B — > 7r form factor 
f + (p 2 ). For brevity, we do not repeat the explicit expressions here. They can be found 
in || together with a detailed discussion of the treatment of ultraviolet and collinear 
divergences. 

The pion distribution amplitude <p n is associated ||, || with the leading twist 2 term 
of the matrix element f\ 

i 

(7r(g) \u{x) llll5 d{0)\0) = -iqj n J du<f^u, p)^ + ... , (13) 

o 

where the ellipses stand for terms of higher twist. The scale dependence of ip n (u, p) follows 
from a Brodsky-Lepage evolution equation. In LO, the solution is given by 

<P*(u,fj.) = ^ s {u) (l + £ a 2n (p)Ct\2u - 1)) , (14) 

where C 2 !^ are Gegenbauer polynomials and a 2n (jj) are multiplicatively renormalizable 
expansion coefficients. While 

<p a w s (u) = Qu(l-u) (15) 

is the asymptotic distribution amplitude, the sum in flT4| ) represents non-asymptotic ef- 
fects with ci2n{n) vanishing at p, — > oo. In NLO flOf , the coefficients a 2n mix under 
renormalization because the Gegenbauer polynomials are no longer eigenfunctions of the 
evolution kernel. A complete NLO expression for <^„.(u, //) can be found in [ |TT|j . 

The spectral density p(si,S2) needed in (|To|) is obtained in LO by substituting the 
hard amplitude (0) and the distribution amplitude ([14]) in @ and taking the imaginary 
parts according to (§). In NLO, one has to add the 0(a s ) hard scattering amplitude T\ as 
specified in (|Tl~l), and use the solution of the NLO evolution equation for <p n . This yields 

p(si, s 2 ) = po(si, s 2 ) + ^— s 2 ) (16) 

47T 



Here, the light-cone gauge, x^A a = 0, is adopted for the gluon field A" 
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with 



1 

Po,i(si,s 2 ) = -U-^lm Sl lm S2 / du ip^u, /i)T 0i i(si, s 2 , u, p) . (17) 



As can be seen from (|T2"|), T only has a single pole in (p + g) 2 . Hence, 

Im S2 T (si,s 2 ,M) = -7iS (l - -^(l -u) - -%u) , (18) 

V m b m b J 

making the convolution with ip n (u,p) very simple. After taking the imaginary part in s\, 
one finds the formal expression for po given in 0. In contrast, T\ also has cuts in (p + q) 2 , 
in addition to the pole. Therefore, Im S2 Ti contains both 5- and ^-functions in s 2 . In the 
following, we use the explicit form given in ||. To obtain the 0(a s ) correction pi(s\, s 2 ) 
one has to convolute Im S2 Ti with the products <^" S C 2 „ . This integration should be carried 
out analytically in order to be able to take the imaginary part in s\. In principle, this 
is possible, by considering the integrals of u k lm S2 Ti for each power k separately. The 
calculation is rather cumbersome, but fortunately not really necessary. It turns out that 
the non-asymptotic effects in p% are much smaller than the overall theoretical uncertainty 
in the LCSR flT0|). It is therefore sufficient to calculate p\ using the asymptotic pion 
distribution amplitude, unless the theoretical accuracy is increased dramatically. 

With the above justification, we single out the asymptotic part of given in flT5| ) 
and perform the integral 

1 

J dw<^ s (w)Im S2 Ti(si, s 2 ,u,n) (19) 


analytically. After taking the imaginary part in si we find 

1 } 

PT( s i' s 2) = -/7r^Im Sl Im S2 y dutp" s (u)T 1 (s 1 ,s 2 ,u,p) 
_ , (r 



2^ + 3ln(-j In (1 + } 



3(3a 3 + 22a 2 + 40a + 24) fa\ ( a\ 3(a 2 + 12a+12) 



a 

-l 22a 2 

2(2 + a) 3 "V2; ' Z V 2 J 4(2 + a) 2 

r ( a 



r \ 1 + r + a 



(20) 




where 



si s 2 st -ml , 01 x 

a=— + — -2, r = 2, (21) 

ml mi s 2 — ml 

X 

and lA 2 (x) = — J yln(l — t) is the Spence function. It is interesting to note that 
contains no In p terms. This is because in T\ the coefficient of In p is given by a convolution 
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of T with the LO kernel of the Brodsky-Lepage evolution equation 0. It thus has to 
vanish when this kernel is folded with the asymptotic distribution amplitude ip® s . 

Finally, the resulting spectral density ( |i~6D is to be inserted in the LCSR fllPl). For 
M\ = Mf = 2M 2 , the leading term coming from p is simply proportional to ip vr (l/2,/i) 
as shown in [ED] . Moreover, the precise form of the duality boundary E is irrelevant in this 
case since only the integration limits at Si = s 2 enter. The NLO correction arising from 
pi, on the other hand, depends on the shape of the integration boundary S. We have 
chosen the triangle s\ + s 2 = 2s^, s$ denoting the threshold of excited and continuum 
states in the 5-channel. This choice allows to carry out the integration over the variable 
r defined in (j2~lD analytically The final LCSR reads 



fBfB'QB'B 7T 



a,C f 



2s s 



M 2 \ e -M*- e -M* ^(1/2, /x) 



47T 



/ / 4- 2 )e^ds + F^\M\mls*^) 

J \ TTli, I 



1m\ 



(22) 



with 



7T , fX 

T + 31n 2 



In 1 + 



x 



3(3x 3 + 22x 2 + 40x + 24) 



In 



x 



2 J 2(2 + x) 3 

3Li 2 (-x) - 3Li 2 (-x - 1) - 3 ln(l + x) ln(2 + x) 
3(3x 2 + 20x + 20) 6x(l + x) ln(l + x) 



+6Li 2 ( -- 



4(2 



x) 



(2 + xf 



(23) 



In ([22]), we have added the contributions F^ 3 '^ from the pion distribution amplitudes of 
twist 3 and 4. The twist 3 contribution is quantitatively as important as the twist 2 term, 
and has to be taken into account in phenomenological applications. Presently, F^ 3 '^ is 
only known in LO JE]. The lack of knowledge of the NLO corrections hinders a complete 
NLO analysis. 

The analogous LCSR for fr>j e D*gD*Dn is obtained from fl22| ) by formally replacing b 
with c, B with D, and B* with D*. 



4. The values of the parameters used in the numerical evaluation of the LCSR (|22|) 
and the two-point sum rules (2SR)f] for the corresponding decay constants are discussed 
in [|H] and ||. A brief recapitulation may therefore suffice here. For the 6-quark mass and 
the corresponding continuum threshold we take 



4.7 ±0.1 GeV, 



35 =f 2 GeV 2 



(24) 



Variations of m& and within the allowed intervals are correlated by the stability re- 
quirement of 2SR as indicated by the opposite signs in (|24]). The renormalization scale 
of a s is put equal to the factorization scale p which will be fixed later. The running 
coupling constant is taken in two-loop approximation with Nf = 4 and =315 MeV 



^In NLO, they are given in 
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corresponding to a s {mz) = 0.118 [13|]. The Borel parameter in LCSR is constrained to 
the interval 6 GeV 2 < M 2 < 12 GeV 2 , where the usual conditions of small (< 10%) twist 
4 corrections and moderate (< 30%) contributions from heavier states are satisfied. In 
2SR the allowed range is 2 GeV 2 < M 2 < 6 GeV 2 . In the charm case, the corresponding 
parameters are given by 

m c = 1.3±0.1 GeV, 5^ = 6^1 GeV 2 , (25) 

and A^L = 380 MeV. The ranges of the Borel parameter are 2 GeV 2 < M 2 < 4 GeV 2 
in LCSR and 1 GeV 2 < M 2 < 2 GeV 2 in 2SR. The meson masses are m B = 5.279 GeV, 
m B . = 5.325 GeV, m D = 1.87 GeV, m D * = 2.01 GeV, and / w = 132 MeV. Finally, for 
the pion distribution amplitude at u — 0.5 we adopt the value ip^(0.5, /i) = 1.2 at ji = 1 
GeV determined in |14] from the LCSR for the pion-nucleon coupling. This value is also 



consistent with the distribution amplitude ip n (u, used in |], || to calculate the form 
factor / + . The values of the vacuum condensates appearing in 2SR can be found in B. 
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Figure 1: Scale dependence of the twist 2 term in the sum rule estimates of (a) fs fis* 9 b* 
and (b) gB*Bn : NLO (solid) in comparison to LO (dashed). 

Let us first discuss the twist 2 term in the LCSR (^) being the main subject of this 
paper. This term is plotted in Fig. la as a function of fi for M 2 = 8 GeV 2 . As can be 
seen, the NLO correction is numerically large, increasing fsfB*9B*B-K by about 50%. It 
has the same sign as the NLO corrections to the 2SR for fs and fg which amount to 
about 30 % and 20%, respectively. Consequently, in the coupling gs*B-K itself, the NLO 
effects cancel almost completely as demonstrated in Fig. lb. A similar cancellation takes 
place in the sum rule calculation of the B — > 7r form factor / + @> 3- As shown in ||, 
the origin of this cancellation can be understood by considering the heavy quark limit of 
the sum rules. In this limit, one has a partial cancellation of large ln(mb/o;o) appearing 
in fsf + and fs, where ojq = (s$ — ml) jlm^ determines the size of the duality interval 
which is mass-independent. In the present case, we observe a complete cancellation of 
these logarithms in the heavy-quark limit of the LCSR ([22|) and the 2SR for fsfs*- As 
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Table 1: Sum rule predictions for B and B* mesons 



LCSR 


f B f B* 9 B* Bit 


0.80±0.23 GeV 2 


2SR 


Sb 


180±30 MeV 


LCSR 
2SR 


fB*9B*Bw 


4.44±0.97 GeV 


2SR 




195±35 MeV 


LCSR 
2SR 


9B*Btt 


22 ± 7 



a result, the NLO correction to the coupling constant 9b*b-k contains no logarithms of 
m b at all. Furthermore, the scale-dependence of f b f b* 9b* b-k is moderate. In LO, it only 
comes from the evolution of the distribution amplitude (p n , while in NLO it also arises 



from the running of a s (fi). Making the usual choice \x = \jm B — ml = 2.4 GeV, a scale 
which is of the order of the Borel mass M, and allowing fi to vary from one half to twice of 
the above value, fsfB*gB*Bn (NLO) varies by 10% relative to the nominal value. Again, 
this is comparable to the NLO scale-dependence of fs and f% so that 9b*b-k is almost 
scale-independent as can be seen in Fig. lb. 

Similar features are observed in the D-meson case. In particular, the NLO correction 
to JdJd*9d*d-k of roughly 15% is compensated by the NLO corrections to jo and fa* 
of about 10% each. It may look surprising that the perturbative corrections are smaller 
in the charm case than in the beauty case, in contrast to the naive expectation from 
asymptotic freedom. The reason becomes again clear in the heavy quark limit: when 
going from beauty to charm, the growth of a s (fi) is overcompensated by the shrinkage of 
the above-mentioned logarithm ln(m;,/ 'ujq). 

The numerical predictions including the twist 3 and 4 contributions in LO are 
summarized in Tables 1 and 2. In dividing out the pseudoscalar and/or vector meson 
decay constants from the LCSR results, one can follow different procedures. The simplest 
possibility is to divide the complete sum rule (^) by the NLO values of these constants. 
Alternatively, one may divide the twist 2 term in (22) by the NLO values of decay con- 
stants, while dividing the twist 3 and 4 terms by the LO values. The numerical difference 
lies within the overall theoretical uncertainty quoted in the Tables. We prefer the first 
option. For 9b*b-k, the uncertainty can be estimated as follows: 

(a) Variation of the Borel parameters in the allowed ranges leads to a variation of the 
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Table 2: Sum rule predictions for D and D* mesons 



LCSR 


f d! d*9d* D-k 


0.54±0.15 GeV 2 


2SR 


fo 


190±20 MeV 


LCSR 
2SR 




2.84±0.55 GeV 


2SR 


Jd* 


270±35 MeV 


LCSR 
2SR 


9D*Dtt 


10.5 ± 3 



coupling by about 10%. 

(b) If the 6-quark mass and the continuum threshold s^ are varied simultaneously 
such that maximum stability of the sum rules for f% and f B is achieved, gB*Bn varies by 
25%. ' ' ' 

(c) The uncertainty coming from the non-asymptotic terms in the pion distribution 
amplitudes amounts at most 10%. This can be checked by discarding the non-asymptotic 
effects altogether. 

(d) For the unknown NLO corrections beyond twist 2 we assign a 15% uncertainty. 
Since twist 3 contributes about 50% to the sum rule fl22| ) this actually assumes an unknown 
correction to twist 3 and beyond as large as 30%. 

(e) The twist 4 contribution to the LCSR (22) is roughly 3%. This can be taken as 
an indication of the uncertainty due to the neglect of twists higher than four. 

(f) The size of the non-asymptotic terms in p implies that they should affect the 
0{a s ) spectral density pi by less than 5%. 

(g) In (|ToD there is some ambiguity in the choice of the duality boundary £ for p\. 
A reasonable way to estimate the maximum uncertainty is to put the upper integration 
limit Sq in (|22|) to infinity which corresponds to no continuum subtraction at all. The 
change is less than 1%. 

(h) The sensitivity to the precise value of the quark condensate density is rather small, 
because of a compensation of the effects in 2SR and in the twist 3 term of LCSR. The 
inaccuracy of higher- dimensional condensates in 2SR is unimportant. 

Adding (a) to (h) in quadrature, the overall uncertainty on gB*B-n is shown in Table 1 
together with the uncertainties in the related quantities estimated in a similar way. The 
corresponding error estimates for the D-meson case are listed in Table 2. 
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A comprehensive comparison with predictions obtained by different methods can be 
found in Table 1 of ref. n. New estimates of gB*Bir in the framework of lattice QCD 



15fl and in a relativistic quark model [16j have appeared very recently. While our result 
shown in Table 1 agrees with the lattice estimate within uncertainties, it is smaller by a 
factor of 2-3 than the quark-model prediction. 

5. As pointed out in the introduction, the B*Btt matrix element determines the 
B — > 7r form factor f + (p 2 ) near the kinematic limit p 2 = (rug — m n ) 2 , where the I?*-pole 
dominates. More definitely, in the single-pole approximation, one has 

/V) = 9 f °: gB 'T, 2 s ( 26 ) 

2m B *{l — p 2 /m z B «) 

with the normalization f B *gB*Bn given in Table 1. Since the NLO corrections to the 
approximation ( [26] ) at large p 2 and to f + (p 2 ) calculated from LCSR directly at low and 
intermediate p 2 are both small, the two descriptions match in NLO as well as they do in 
LO [[IJ]. This also holds for the D-meson form factor. 

Furthermore, from go*D-K given in Table 2, one obtains : 

T(D*+ -> D°ir + ) = I q I 3 = 23 ± 13 keV, (27) 



q being the decay 3-momentum. The current experimental limit ]13| 



r(/r+ -> D°n + ) < 89 keV (28) 

is still too high to challenge the theoretical prediction. 

In conclusion, we have made some progress in the NLO analysis of the LCSR for the 
B*Bn and D*Dtt coupling constants. The perturbative 0(a s )-correction to the twist 2 
term of this sum rule has been calculated removing one of the main uncertainties. The 
remaining uncertainties in gs*B-K and gD*Dn are mainly due to the heavy quark masses, 
the uncertainty in the pion distribution amplitudes of twist 2 and 3 at u ~ 1/2, and the 
lack of the NLO corrections to the twist 3 contributions. Thus, there is still plenty of 
room for improvement. 
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